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Abstract 

A one dimensional array of x-ray detection pixels is formed 
by orienting a silicon strip detector edge on, with its strips 
parallel to the x-ray beam. Spatial resolution is governed 
by the pixel size, defined by the silicon thickness and strip 
pitch The detector orientation presents a large apparent 
depth for. x-ray conversion, resulting in a high quantum 
efficiency. Coarse, but very stable x-ray energy determina- 
tion is made by counting the number of conversions as a 
function of depth in each pixel: The strips (pixels) are sub- 
divided into segments, each sampling a different conversion 
interval. Custom CMOS ICs provides each detector ele- 
ment with an independent low noise preamp, shaper, com- 
parator and 16-bit counter. During the acquisition phase, 
each channel counts x-ray conversions, at rates in the 1 
MHz range. Readout consists of uploading the counter 
sums (8 bytes/pixel) to the host computer. 

I. INTRODUCTION 

The imaging task in transmission radiography is to mea- 
sure the differential x-ray attenuation of a subject across 
its area. Conventional. film based radiography and many 
of the so called "digital" techniques are based on energy 
conversion (e.g. x-rays to visible light photons) followed 
by an integration stage. These methods produce a signal 
proportional to the total energy of the transmitted beam 
deposited in the converter. For a monoenergetic beam, this 
is proportional to the number of transmitted x-rays. For a 
polyenergetic beam, the signal is still proportional to the 
deposited energy of the transmitted beam. However, this 
is no longer proportional to the number of transmitted 
x-rays, hence hence does not represent a direct measure 
of the attenuation. This effect and the stochastic process 
governing the energy conversion [1], reduce the statistical 
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significance of measured differential signals. A high effi- 
ciency, direct quantum detection of individual x-rays does 
not suffer from either of these effects and is capable of 
producing images with correspondingly higher statistical 
significance for a given patient dose. In this paper, we 
describe such a system under development at LBNL. 

Fully depleted, high resistivity segmented (strip or pixel) 
silicon diode arrays are in widespread use as position sen- 
sitive detectors for charged particles in elementary particle 
and nuclear physics experiments. Large areas are instru- 
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Figure 1: edge on detector under development 

mented by orienting planar arrays normal to the particle 
direction. While 300 microns of silicon is adequate for de- 
tection of minimum ionising particles, it is not well suited 
for x-ray detection above 10-15 keV (see Table 1.) How- 
ever, by orienting the detector array "edge on" (Figure 1) 



the x-rays see much more material, offering a much higher 
detection efficiency. An edge on strip detector becomes 
a 1-dimensional pixel array with pixel size defined be the 
strip pitch and detector thickness. A 2-dimensional image 
is then formed by scanning the subject with a fan beam. 
The fan beam configuration results in enormous reduction 
in scatter without the efficiency loss associated with an 
anti-scatter grids between the subject and detector. 

The concept of orienting a planar detector edge on in 
order to overcome the efficiency problem is not new[2, 3]. 
However, our concept turns the poor attenuation proper- 
ties of silicon to an advantage: The pixels are segmented 
in depth (Figure 1) and the number of x-rays interacting 
in each segment is counted. The measured x-ray penetra- 
tion profile provides energy information of the emerging 
beam. Table 1 lists x-ray attenuation lengths in silicon for 
several energies used in the diagnostic radiology [4]. These 

Table 1: X-ray attenuation lengths in silicon 



E (keV) 


20 


30 


40 


60 


80 


100 


/i" 1 (mm) 


0.97 


3.01 


6.15 


13.40 


19.17 


23.24 



relatively large distances allow the readout electronics to 
be spaced over corresponding distances, simplifying sys- 
tem design. The resulting energy information is robust; 
it depends only on the underlying attenuation properties 
of x-rays in silicon, not channel-to-channel variations in 
leakage currents, gains, and energy thresholds. The result- 
ing energy information can be used to manipulate tissue 
contrast difference in the same way that dual energy tech- 
niques do, without any image registration problems[5]. 



II. DETECTOR LAYOUT 

The optimal pitch, length and segmentation of the detec- 
tor strips will depend on the imaging task. If total number 
of strips N > L/D, where L is the x-ray source- detector 
distance and D is the strip length, the effect of parallax 
requires that the strips point back to the source. Thick 
subjects requiring higher energy beams will require longer 
strips to remain efficient, relative to thin subjects. Very 
high hit rates can be accommodated by simply increasing 
the strip segmentation. The detectors designed and being 
fabricated at LBNL are not optimized for any particular 
imaging task. They are 300 /im thick n-type high resistiv- 
ity silicon and have 128 45 mm long strips. Each is divided 
into four segments of length 3, 6, 12 and 24 mm(upstream 
to downstream.) The strips have a 100 /im pitch at the 
upstream end of the detector and point toward the x-ray 
source 1.15 m away. A set of collimators upstream from 
the subject define the fan beam; a downstream collimator 
limits detector acceptance to transmitted, not scattered 
radiation. 
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Figure 2: Front-end architecture 

III. FRONT-END 
ELECTRONICS 



Figure 2 illustrates the architecture of the front-end read- 
out chip which is being developed using the HP 0.8/im 3 
metal process available through MOSIS. Each channel of 
the analog portion consists of an integrator and preamp 
with a 100 ns peaking time and gain of 860 mV/fC. A fast 
preamp reset allows operation at 1 MHz with an equivalent 
input noise of 200 r" for detector capacitance Co ~ 1 - 3 
pf. Each preamp is DC coupled to the input of a compara- 
tor. The RMS spread in shaper outputs levels measured 
on 17 64-channel chips was 8 mV, corresponding to 58 e". 
Setting the discriminator at 5 a should result in noiseless 
operation with an energy threshold of 3.6 keV. 

The digital section consists of one 16-bit counter per 
channel. At the start of a data acquisition cycle, a RESET 
is issued which clears the counters. At the end of the ac- 
quisition phase, a PARALLEL LOAD is issued causing the 
contents of the counters to be transferred into a 1024 bit 
(64 channels times 16 bits) shift register. The 64 channels 
of digital data are readout serially, taking about 50 /is. In 
principle, once the counter contents have been transferred 
to the shift register, a new acquisition phase could begin, 
keeping the deadtime at a few hundred ns/readout cycle. 
The shift registers have a DATA IN line for daisy chaining 
multiple chips. together. 
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IV. EXPECTED 
PERFORMANCE 

At the time our detectors were being designed, we had not 
attempted to optimize for a specific imaging task; our goal 
has been to provide a system to demonstrate proof of prin- 
ciple. It is useful to simulate the behavior of the system 
for a concrete, demanding imaging task. We chose mam- 
mography as an application since the compressed breast 
is a relatively simple and well defined subject to model. 
Other imaging tasks may be well suited, but more difficult 
to simulate. For example, chest imaging could certainly 
benefit from the linear, virtually unlimited dynamic range 
. of an x-ray counting detector; bone density measurements 
could utilize the energy sensitivity for soft tissue cancella- 
tion, etc. 

A. Dose Calculation 

In any radiographic procedure, an assessment of the risk 
that the radiation poses to the patient must be made. In 
mammography, the mean glandular dose (MGD) is the pre- 
ferred quantity used to assess risk[6]. We would like to 
express the image quality, i.e. the statistical significance 
of differential attenuation measurements, as a function of 
patient dose, which depends on incident x-ray energy spec- 
trum, breast composition and thickness. A detailed x- 
ray transport Monte Carlo was developed [4] to calculate 
Dc(E\b) y the mean glandular dose per incident x-ray as a 
function of energy (E) and breast composition/thickness 
(6). For a given breast (6) and incident x-ray spectrum, 
*(£), MGD=fdE&(E)D c (E t b). 

Since the fan beam is expected to scan the entire surface 
area of the breast, x-rays at each energy were distributed 
randomly over the entire area. The breast is modelled 
as a set of slabs, each a uniform mixture of that mate- 
rial's atomic constituents, stacked parallel to the beam. 
Two 0.4 cm thick adipose layers representing skin surround 
a central mixture of glandular (50%) and adipose (50%) 
tissues[7]. Several total thickness were simulated ranging 
from 4 to 9 cm, but here results are restricted to 5 cm. At 
the top and bottom of the stack are 0.3 cm Lexan compres- 
sion paddles. The transverse size of the stack was a rect- 
angle 12.7 x 25.4 cm 2 . 10 6 x-rays were generated for each 
energy; each was propagated through the different volumes 
and allowed to undergo coherent, incoherent and/or photo- 
electric interactions. Each x-ray was transported through 
the material until it escaped or underwent. a photoelectric 
interaction. D G (E) was calculated by dividing the total 
energy deposited in the central breast region by the mass 
of that region and dividing the result by 10 6 . 

The Monte Carlo also yielded the scatter/primary ratio, 
S pt as a function of breast thickness, x-ray energy and 
collimation. For an uncolHmated beam, S p varied between 
0.5 and 1.1 as the breast thickness varied between 4 and 
9 cm, showing little energy dependance. For the fan beam, 
S p was more then a 100 times smaller. 
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B. Efficency 

Figure 3 shows the calculated quantum detection efficiency 
for a 45 mm long detector. The dashed curve shows the 
total interaction probability while the solid curves show 
the interaction probabilities for inelastic interactions which 
deposit energy above the listed thresholds. The quantum 
efficiency of the system will exceed 67% over the entire 
diagnostic range for a discriminator threshold of 3.6 keV. 
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Figure 3: Quantum detection efficiency 
C. Image Signal/Noise 

The signal to noise ratio (S/N) for a fixed tissue dose 
is often employed as a figure of merit in the evaluation 
of an imaging system. Consider a uniform target region 
imbedded within a uniform background. In this case, 
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Figure 4: Detectibility of a microcalcification in mammog- 
raphy; monoenergetic beam 

S/N= \N t -N b \/y/N t + N ht where N t and N b are the mea- 
sured numbers of transmitted x-rays in the target area and 
background region of equal area. In general, N t and in- 
clude scattered as well as primary transmitted x-rays. For 



a subject scatter to primary ratio, S p> the statistical sig- 
nificance of a different ial measurement is decreased by ap- 
proximately yjl + S P for a. given patient dose. Our x-ray 
transport Monte Carlo shows that the fan beam configu- 
ration results in S p < .01 and we ignore it. 

We have evaluated the the S/N factor for a mammog- 
raphy task: the detection of a microcalcification in a uni- 
form background field. The calcification is taken to be 30 
mg/cm 2 of hydroxyapatite [Ca 5 (P0 4 ) 3 OH] displacing the 
central 250 ^m of a 5 cm thick compressed breast. The 
solid curve in Figure 4 shows the S/N for our detector 
(threshold = 3.6 keV) as a function of x-ray energy for 
a fixed mean glandular dose of 1 mGy and target size of 
3 x 10" 4 cm. (for a given spectral shape, S/N scales with 
the the square root of the dose x area.) The dashed curve 
corresponds to a fan beam imaging system employing a 
typical mammographic intensifying screen (31.7 mg/cm 2 
of Gd 2 0 2 S) perfectly coupled to an ideal visible light inte- 
grator [8]. 

Figure 5 shows, for the same task, the S/N comparison 
for the two imaging methods as a function of peak kilovolt- 
age settings of a simulated mammographic x-ray source 
(Mo anode, 10° target angle, 0.8 mm Be window plus 
added filter. The generated spectra [9] were normalized to 
1 mGy mean glandular dose) The solid curves correspond 
to 30 /im of added Mo filtration. The dotted, dashed and 
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Figure 5: Detectibility of a microcalcification in mammog- 
raphy; Mo anode spectrum 

dot-dashed curves correspond to 0.25, 0.50 and 1.00 mm 
added Al in place of the Mo filtration. It should be noted 
that our method has a higher S/N and a shallower fall 
off with increasing kilovoltage. The Al filtered spectra are 
harder than the Mo filtered spectrum with the thicker fil- 
ters showing superior S/N, but require higher mAs values. 
Our detector achieves the same optimal S/N for the Mo 
or 0.25 mm Al filtered spectra, though the later achieved 
it at 60% of the tube loading. The peak position in Fig- 
ure 4 suggests that a Rh target (fluorescent lines at 20 and 
23 keV as opposed to the Mo lines at 17 and 209 keV) 
would be better suited for our detector. 



D. Energy Content 

Figure 6 shows the simulated energy spectrum of detected 
x-rays in the first strip (a) and the third + fourth strips 
(b) resulting from an unattenuated, W anode, Al filtered, 
120 KV p x-ray source. The mean energies of the two spec- 




Figure 6: Energy spectrum of detected events in (a) strip 
1 and (b) strips 3+4 

tra are 34 keV 68 keV respectively, illustrating the energy 
sensitivity of the system. We anticipate the ability to ma- 
nipulate image contrast using dual energy basis decompo- 
sition technique[5]. As an example, image clutter resulting 
from the variation in thicknesses of two tissue types can 
interfere with the task of measuring a third, overlapping 
tissue type. This clutter can be removed by forming a 
linear combinations of images acquired at two different en- 
ergies: T = Ti cos</> LA + T H s\n<f> LAy where 7} and Tfi are, 
for each pixel, the log attenuations measured at low and 
high energy, T is the resulting image projection, and <f> LA 
is the material look-alike projection angle of reference^)]. 
In general, this projection decreases the S/N for the third 
tissue. We have simulated several different x-ray spectra 
(Mo and W targets, various filters) and strip segmenta- 
tions to evaluate the detectibility of microcalcifications in 
breast tissue after removing clutter produced by adipose 
and glandular tissue. Our initial conclusion is that, for this 
application, the loss in S/N is too severe to offset the gains 
of clutter removal. We are presently pursuing other pos- 
sible applications, such as chest imaging or bone density 
measurements. 

V. Conclusion 

We have presented a detection scheme for counting indi- 
vidual x-rays with high quantum efficiency using a silicon 
strip detector presented edge on to a fan beam. Robust 
energy sensitivity is achieved by measuring the depth de- 
pendance of the x-ray conversions. Very high system band- 
width is achieved by counting and storing the x-ray inter- 
action sums on the front end electronics. The detectibility 
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of microcalcifications in breast imaging, expressed as the 
image S/N is expected to surpass that of screen based sys- 
tems. We are exploring other imaging tasks that could 
benefit from the dual energy capability. 
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